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ABSTRACT 

We describe a program to identify optical counterparts to radio sources from the VLA FIRST survey 
using the Cambridge APM scans of the POSS-I plates. We use radio observations covering 4150 deg^ of 
the north Calactic cap to a 20 cm flux density threshold of 1.0 mJy; the 382,892 sources detected all have 
positional uncertainties of < 1" (radius of 90% confidence). Our description of the APM catalog, derived 
from the 148 POSS-I O and E plates covering this region, includes an assessment of its astrometric and 
photometric accuracy, a photometric recalibration using the Minnesota APS catalog, a discussion of 
the classification algorithm, and quantitative tests of the catalog's reliability and completeness. We go 
on to show how the use of FIRST sources as astrometric standards allows us to improve the absolute 
astrometry of the POSS plates by nearly an order of magnitude to ~ 0.15" rms. Matching the radio and 
optical catalogs yields counterparts for over 70,000 radio sources; we include detailed discussions of the 
reliability and completeness of these identifications as a function of optical and radio morphology, optical 
magnitude and color, and radio flux density. An analysis of the problem of radio sources with complex 
morphologies {e.g., double-lobed radio galaxies) is included. We conclude with a brief discussion of the 
source classes represented among the radio sources with identified counterparts. 

Subject headings: surveys - catalogs - quasars:general - galaxies: general - radio continuum:general 



1. INTRODUCTION 

The optical identification of the first discrete extrater- 
restrial radio source occurred as a result of a telephone 
call from J. S. Hey to the Royal Greenwich Observatory 
on the afternoon of February 28, 1942. Recognizing that 
the source of extensive jamming of British radar over the 
previous two days appeared to follow the Sun, Hey was 
delighted to learn that an unusually large sunspot had 
just transited the solar disk; despite its skeptical reception 
by his superiors. Key's identification proved correct (Hey 
1973). 

In the ensuing decade, progress in the detection of new 
extrasolar radio emitters far outstripped the ability of as- 
tronomers to associate them with optical counterparts. 
The first breakthrough came in 1949 when Bolton, Stanley 
and Slee (1949) identified the Crab Nebula, M87, and NCC 
5128 (Cen A) with three of the brightest radio sources in 
the sky, although they concluded that the bizarre morphol- 
ogy of the latter generally favored a Galactic interpretation 
for radio emitters since "the probability of [such] an un- 
usual object in our own Galaxy seems greater than a large 
accumulation of such objects at a great distance." The fol- 
lowing year, Ryle, Smith, and Elsmore (1950) concurred 



in this conclusion despite finding 0/146 bright {V <4.0) 
stars, 0/21 novae, 0/38 planetary nebulae, 0/29 diffuse 
Galactic nebulae, and 4/5 of the brightest galaxies coin- 
cident with entries in their fifty-source radio catalog. It 
was not until the classic papers of Baade & Minkowski 
(1954a, b), which among other things pronounced Cygnus 
A "an extragalactic affair" , that the era of extragalactic 
radio source identification can be said to have begun. 

The largest radio catalogs in existence prior to 1995 con- 
tained, in total, approximately 100,000 distinct entries. In 
striking contrast to the earliest speculations, fewer than 
20 of these relatively bright radio sources are identified as 
stars. Indeed, fewer than 1000 stellar radio detections have 
been made despite decades of sensitive, targeted searches 
(Hjellming 1988 and references therein; Wendker 1995), 
and < 5% of all cataloged radio sources are Galactic ob- 
jects. A search of the NED database, however, suggests 
that the fraction of identified extragalactic radio emitters 
today is little better than it was in 1950, when 7 out of 
67 known radio sources had identified counterparts (Baade 
and Minkowski 1954b). The problem now is the same as 
it was 50 years ago: the angular resolution of large-area 
radio surveys is generally too poor (several arcminutes) 
to allow for the unambiguous association of cataloged ob- 
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jccts with individual optical counterparts which, at the 
hmit of the POSS-I plates, number >4000 deg~^ at high 
Galactic latitudes. Interferometric surveys at centimeter 
wavelengths can achieve the desired positional accuracy of 
~ 1" but have, until recently, covered only ~ 50 deg^ of 
sky, resulting in fewer than 500 optical identifications for 
radio-selected objects at faint flux levels (see Table 5). 

While considerable information on the class and emis- 
sion mechanism of a radio source can be derived from ob- 
servations of its radio morphology, its spectrum, and its 
polarization characteristics, optical observations are still 
required to establish the source's distance and to classify 
it unambiguously. Eight years ago, we began to construct 
Faint Images of the Radio Sky at Twenty-cm (FIRST) 
with the primary goal of obtaining a very large sample of 
radio sources with positions sufficiently accurate that the 
majority of objects detected could be easily identified on 
the basis of positional coincidence alone. FIRST has been 
designed to cover the same 10,000 deg^ region as the Sloan 
Digital Sky Survey (SDSS), which will obtain deep optical 
images of the northern sky in five colors, and will take spec- 
tra of a million objects over the next decade (Gunn and 
Knapp 1993). However, an archive of the optical sky in 
the form of the National Geographic-Palomar Observatory 
Sky Survey (POSS) plates already exists. The identifica- 
tion of even ~20% of all FIRST sources with counterparts 
at or above the POSS plate limit will immediately provide 
samples of various radio source populations from one to 
three orders of magnitude larger than those in existence, 
and will advance significantly our knowledge of the radio 
universe. 

In this paper, we present the results of an optical iden- 
tification program for FIRST radio sources based on the 
Cambridge Automated Plate Measuring Machine (APM) 
scans of the POSS I plates (McMahon and Irwin 1992). In 
§2, we describe briefly the parameters of the FIRSTsuiyey 
and the catalog derived therefrom. In §3, we describe the 
genesis of the APM catalog and discuss its astrometric, 
photometric, and source classification uncertainties with 
particular attention to those which either 1) the FIRST re- 
sults can help to refine, or 2) arc particularly relevant to ra- 
dio source identification; the following section (§4) demon- 
strates the utility of FIRST as an astrometric calibrator 
by deriving both intraplate and plate-to-plate corrections 
to the optical astrometric solution. Section 5 describes 
the matching algorithms employed and discusses such is- 
sues as false match rates, reliability and completeness of 
the proposed identifications, and the effects of optical and 
radio source morphology. Section 6 discusses color, magni- 
tude, and radio flux density distributions for the more than 
70,000 counterparts we have identifled, while the following 
section (§7) describes the format of the counterpart cata- 
log (the full contents of which are available on the Web). 
In the discussion (§8), we place in context this identifica- 
tion program for faint radio sources, and summarize the 
statistics of the populations represented. A reprise of our 
results and a precis of future work concludes our report 
(§9)- 



2. THE FIRST SURVEY 

In its initial five years, the FIRST project surveyed 
~4670 deg^ of the high Galactic latitude sky to a point- 
source flux density limit of 1.0 mJy at a wavelength of 
20 cm. Using the B-configuration of the Very Large Array 
(VLA^) in bandwidth-synthesis mode, we have produced 
nearly uniform-sensitivity images of the sky with a median 
rms of 0.14 mJy and an angular resolution of 5"; the po- 
sitions for each of the half million sources in our current 
catalog are accurate to better than 1" (90% confidence). 
Details of the observing strategy and the pipeline process- 
ing system that creates the final images from the raw uv 
data are presented in Becker, White, and Helfand (1995; 
hereafter, BWH). 

For the program described here, we have utilized the 
98Fcb04 version of the FIRST catalog, which is available 
from our Web homepage (http://sundog.stsci.edu). The 
catalog contains positions, peak and integrated flux den- 
sities, source morphological parameters, and information 
on the fleld from which each entry was derived. A detailed 
description of the catalog's construction can be foimd in 
White et al. (1997; hereafter WBHG); an exhaustive dis- 
cussion of the astrometric, photometric, and morphologi- 
cal uncertainties is included in that paper and in BWH, 
and is summarized on our Web homepage. 

The catalog has a mean surface density of ^90 "som-ces" 
per square degree. However, it is important to note that a 
signiflcant fraction of radio sources at these flux densities 
are extended. In particular, classical doublc-lobcd radio 
sources in which the components are separated by more 
than ~3" are typically represented by two catalog entries; 
for complex sources extended on scales of ~1', several com- 
ponents can be required to fit adequately the source sur- 
face brightness distribution (see WBHG). We have not yet 
completed a definitive study of multi-component sources; 
even when this is done, it will not be possible to decide un- 
ambiguously whether or not two nearby catalog entries are 
parts of the same source hosted by a single optical counter- 
part, or whether they are chance alignments of unrelated 
objects. We briefiy investigate these issues in §5.4; prior to 
that discussion, we use the term "sources" synonymously 
with "radio catalog entries" . 

For clarity, we present in this paper only the results of 
the APM identification program for the ~ 4150 deg^ in the 
north Galactic cap. This region has approximate bound- 
aries of -1-22° < ^ < -1-58° and 07'^30™ < RA < 17'^30™ 
and includes ^ 382,892 radio sources. Subsequent ex- 
tensions to this work for the two southern cap strips 
(+1.6° > 6 < -2.5° and -7° > S > -11° in the range 
2l'^30'" < RA< OS'^OO'^ for a total of 610 deg^) as weU 
as the remainder of the north cap region will be added to 
the Web site as they become available. 



3. THE APM CATALOG 

3.1. Introduction 



It has been almost 50 years since the National Geo- 
graphic Society-Palomar Observatory Sky Survey (POSS) 

^ The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooperative agreement with the National 
Science Foundation. 
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was carried out (1949 to 1958) using the 48" Oschin 
Schmidt at Palomar Observatory (Harrington 1952; 
Minkowski and Abell 1963; Lund and Dixon 1973). Since 
a second-epoch Palomar Survey is underway, the original 
POSS is now usually referred to as the first-epoch POSS or 
POSS-I to differentiate it from POSS-II (Reid et al. 1991). 
The survey covers the 31,000 deg^ of sky north of 5 ~ —33° 
with 6.4x6.4 degree plates in two colors: blue (O), cen- 
tered at ~ 4100 A covering the range between 4900 A and 
the atmospheric cutoff at ^^3300, and red (E), a 600 A 
wide band centered at ^^6500 A. Plots of the POSS-I band- 
passes compared with the standard Johnson photometric 
system can be seen in Evans (1989). The nearest Johnson 
system bands to O and E are B and R respectively. We 
provide color equations for these bands in section 3.3. 

Between 1991 and 1995, the SERC Automated Plate 
Measuring Alachinc (APM) at the Institute for Astron- 
omy, Cambridge (Kibblcwhitc ct al. 1984) was used to 
digitize these plates at a resolution of 7.5 microns per pixel 
(0.5"), the highest spatial resolution yet applied to these 
images (McMahon and Irwin 1992). An object catalog 
has been constructed from these data which includes all 
objects down to the plate limits — ~20.0 in E and ~21.5 
in O — and contains approximately 2000 stars and 2000 
galaxies deg^^ at high Galactic latitudes. The catalog con- 
tains positions, magnitudes, morphological classification 
parameters, major and minor axes, and isophotal areas 
for each source; a merged catalog which matches objects 
between plates also contains a color (or an upper limit 
thereto) for each entry. An automated classification algo- 
rithm interprets the morphological parameters to classify 
each object. Here, we present the basic procedures which 
establish the astrometric and photometric calibration of 
the APM catalog, and discuss the limits of the image clas- 
sification system. 

3.2. APM Astrometry 

The APM machine measures the x and y positions of 
all objects detected. The conversion relationship between 

these measured APM positions and celestial coordinates is 
derived by matching stars in the Tycho-ACT catalog (Hog 
et al. 1997; Urban, Corbin & Wycoff 1998) with stars de- 
tected on each plate using a 'standard' six plate-constant 
model that allows for shift, rotation, scale, and shear. The 
algorithm uses iterative 2.5(t clipping to give a robust fit; 
the typical rms on the fitted positions of the Tycho stars 
are 0.4"-0.8". 

Irwin (1994) has studied the two-dimensional systematic 
errors in an earlier version of the APM catalog positions 
by investigating the intraplatc residuals between the mea- 
sured positions for bright stars in the Positions and Proper 
Motions Catalog (Roser and Bastian 1991) and the astro- 
metric fit. He found significant, systematic residuals rang- 
ing up to 0.5". In the version of the APM catalogue used 
herein, the astrometric analysis uses the more recent ACT 
catalog (Urban, Corbin & Wycoff 1998). 

A residual map generated from this analysis is applied 
to positions in the standard APM catalog available at 
http://www.ast. cam. ac.uk/^apmcat. We discuss below 
(§4) the significant improvement in overall astrometric ac- 
curacy that can be derived from a comparison of radio 
source counterparts in the APM with FIRST smyey source 



positions. 

3.3. APM Photometry 

The APM measures photographic density rather than 
flux; moreover, the central regions of all objects more than 
a factor of ~10 brighter than the sky produce a nonlinear 
response and/or are saturated. The algorithms used to 
overcome these inherent difficulties are discussed in detail 
by Irwin (1985). Briefly, a local background is determined 
for each of ^500000 locations on each plate by producing 
a histogram of the pixel values in 64 x 64-pixel regions 
(32" X 32") and finding the mode of each distribution; 
two-dimensional smoothing is applied to these half mil- 
lion background estimates to derive a background model 
for the plate. The image detection algorithm then finds 
connected regions of pixels above a threshold level (typi- 
cally 2(7 above the estimated background level for the given 
plate position). This background- following technique has 
the advantage that faint objects lying in the halos of bright 
objects can be detected. However, large objects such as 
bright stars and galaxies with angular extents > 30" have 
their raw fluxes underestimated. An additional problem 
for large images is that the limited memory available to 
the software means that bright objects sometimes over- 
flow the pixel buffers and are lost. This occurs for images 
with sizes greater than roughly 1-2 mm (i.e., 1-2'), corre- 
sponding to stellar magnitudes brighter than ^ 9.0. 

Another inherent problem arises in attempting to derive 
magnitude estimates for extended objects from saturated 
images. Saturation effects can be corrected for in stars by 
assuming that stellar images have an intrinsic density pro- 
flle independent of magnitude, and that this profile can be 
derived from the unsaturated parts of stellar profiles. A 
high signal-to-noisc intrinsic profile is constructed by tak- 
ing the core from faint stars and the wings from brighter 
stars (see Bunclark & Irwin (1983) for further details). 
This profile can then be integrated and used to derive a 
calibration curve to convert saturated stellar magnitudes 
to a linear system. In the APM catalog, this calibration 
is applied to all images. This has the unfortunate conse- 
quence that galaxies, which have shallower surface bright- 
ness profiles and lower central surface brightnesses than 
stars of the same total magnitude, will have their magni- 
tudes over-corrected. This is a fundamental problem for 
galaxy photometry determined from photographic sky sur- 
vey plates (sec Metcalfe, Fong, and Shanks (1995) for a 
discussion) . 

The basic APM catalog is defined to have a red-band 
(E) plate limit of to(Rc) = 20.0. This limit was estab- 
lished during the early stages (^1991) of the creation of 
the APM catalog via comparison with ^10 photometric 
sequences (Evans 1989; Humphreys et al. 1991). Simi- 
larly, a single slope of 1.10, was assumed in converting be- 
tween the linearized APM magnitudes (Bunclark & Irwin 
1983) and the a Lyrae-based Johnson magnitude system. 
It was noted at the time that there were significant devia- 
tions (~1 mag) from a simple linear relation at magnitudes 
brighter than ~15. This is not surprising bearing in mind 
that the POSS-I glass plates measured by the APM are 
copies that may have different degrees of saturation and 
have had their contrast stretched to enhance faint features. 
The assumption of a constant flux limit seemed reasonable, 
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since the plates were all taken in similar dark sky observ- 
ing conditions with exposure times that were adjusted to 
ensure uniform sensitivity. A similar assumption is made 
in all modern photographic cameras where it is assumed 
that all photographic film has the specified speed. The 
blue band (O) limit was defined with respect to the red 
limit; for the 428 fields available in March 1999, this has 
a range of m(Bo)=20.6-21.3 (ilcr). 

Eventually, a full photometric recalibration of the APM 
using the Guide Star Photometric Catalog (GSPC Post- 
man et al. 1998a) CCD sequences is planned. Preliminary 
comparisons with CCD photometry for ~ 5% of the POSS- 
I plates show that the APM magnitudes for stellar objects 
have a global rms uncertainty of 0.5 magnitudes over the 
range 16 to 20, the range in which most FIRST coun- 
terparts lie. As discussed above, the uncertainties in the 
magnitudes of galaxies are more complex, since galaxies 
have a range of surface brightness distributions, and hence 
may have complex, partially saturated surface brightness 
profiles on the POSS-I plates. This is compounded by 
the range in calibration slopes observed. At faint magni- 
tudes (18-20) where the image profiles are unsaturated, 
the APM magnitudes may be more reliable, but it is left 
to the reader to verify this where precise magnitudes are 
required. For many programs, a uniform set of magnitudes 
or uniform selection criteria are more critical. 

It is also worth noting that almost 50 years has elapsed 
between the epochs of the POSS-I and FIRST surveys, so 
that optical variability is an additional uncertainty. Hook 
et al. (1994) have studied the long term variability of ra- 
dio quiet quasars and found that over a rest-frame period 
of ~10 years, a typical m{B) = 19 quasar varies by 0.20 
magnitudes (rms). The longer-term variability of FIRST 
optical counterparts could be studied via a comparison be- 
tween the POSS-I plates and the POSS-II or UKST plates. 
A CCD investigation of the long-term variability of ^ 200 
quasars from the FIRST Bright Quasar Survey (White et 
al. 2000) is reported elsewhere (Helfand et al. 2001). 

It is useful to be able to convert the magnitudes of the 
O and E bands to the nearest bands in the Johnson Vega- 
normalised magnitude system. Evans (1989) has found 
that over the color range 0.0 < B— R < 1.5, 

R- E = {om±om){B - R) (1) 

Thus one can assume; 

R = E (2) 
Evans (1989) also found 

0-£;= (1.135 ±0.035)(B-i?) (3) 
Thus we have; 

B-R = 0.88(0 - E) (4) 
Finally it follows from (1) and (3) that: 

B = O-0.12{O-E) (5) 

Assuming central wavelengths of 4100A and 6500A and 
the Hayes and Latham (1975) calibration of Vega the 

^ The APS databases are supported by the National Science Foundation, 
of Minnesota, and are available at http://aps.umn.edu/. 



monochromatic zero point of the O and E bands, i.e., the 
flux corresponding to a magnitude of zero, is 4550 Jy and 
2980 Jy respectively. Caution is advised when using these 
conversions since the presence of emission lines in the O 
and E filters will eflFect the conversion. 



3.4. Photometric Calibration Using APS 

In an attempt to improve the APM photometric accu- 
racy and uniformity, the APM magnitudes in regions cov- 
ered by the FIRST smvey were recalibrated plate- by-plate 
using magnitudes from the Minnesota Automated Plate 
Scanner POSS-I catalog (APS^, Pennington et al. 1993) 
which are more uniform than the APM magnitudes be- 
cause they were calibrated on a plate-by-plate basis. The 
APS catalog was created by scanning the same POSS-I 
plate material as the APM and so should be fully con- 
sistent with the APM since it has the same bandpasses, 
epoch, and so on. 

An alternative approach would have been simply to use 
the APS catalog in place of the APM catalog to get opti- 
cal counterparts for the FIRST sources. We preferred the 
APM catalog because the APS catalog is not complete 
over the POSS-I area and does not cover the southern sky 
at all, and because the APS catalog retains only sources 
that appear on both the red and blue plates, discarding a 
significant fraction (-^40%) of the faint radio source coun- 
terparts near the plate limits. Analysis of the reliability 
of these single band detections is discussed in section 5.2 
and shows that 92.7% of the blue-only matches within 1" 
are real, and 97.5% of the red-only matches within 1" are 
real matches. 

We matched the entire FIRST catalog against the APS 
catalog, extracting all objects within 20" of each radio 
position, and then matched the resulting list of optical 
sources with the equivalent APM/FIHST match fist. For 
each APM plate, we determined independent linear fits 
(zero-point and slope) for E and O that transform the 
APM magnitudes to the APS scale. The fit was based 
on APM/ APS matches closer than 5" that are classified 
as stellar by the APS (which uses a different photometry 
method for non-stellar objects.) 

For 8 of the 148 POS'S-l plates covering the FIRST 
area, there were insufficient APS sources available to de- 
termine the photometric calibration because the corre- 
sponding plate was unavailable in the APS catalog. For 
those plates, we bootstrapped a photometric solution us- 
ing APM sources in the overlapping regions of neighboring 
plates. The set of zero-points and slopes for all 148 plates 
is available on the FIRST website. 

This calibration procedure substantially improves the 
APM photometry. This is clearly seen in the magnitude 
discrepancies for APM objects in the plate overlap regions, 
which provide two or more independent magnitude mea- 
surements per object. Figure 1 shows the distribution of 
APM magnitude differences before and after calibration; 
the rms scatter decreases from 0.45 to 0.30 magnitudes, 
and the scatter for bright sources is reduced by an even 
larger factor. Since the plate overlap regions lie at the ex- 
treme edges of the POSS-I plates, they are probably the 

the National Aeronautics and Space Administration, and the University 



Optical IDs for 70,000 Radio Sources 



5 



worst-calibrated areas on the plates; consequently, we esti- 
mate that the recalibrated APM magnitudes are accurate 
to better than ~ 0.2 magnitudes rms. 

Figure 2 shows that the principal problem with the APM 
magnitudes is a magnitude-dependent error. We display 
the differences between the mean calibrated O and E mag- 
nitudes and the original APM magnitudes as a function of 
magnitude. Both colors display a quasi-linear trend, with 
the total error varying by ~ 0.5 magnitudes over an eight- 
magnitude span. While these curves have been derived 
from a small subset of all APM scans, it is likely that 
application of the corrections they imply will improve sig- 
nificantly the photometric accuracy of the catalog. 

3.5. APM Classification 

The APM scans result in a parameterization of each 
detected image which includes an x,y position, a peak 
intensity, a total isophotal intensity, second moments of 

the intensity distribution, and arcal profiles (defined as 
the number of pixels above preset levels which increase 
by powers of two above the threshold level). In addition, 
a classification parameter is calculated which reports by 
how many sigma the object differs from the stellar point- 
spread function on each plate. The stellar psf is derived 
as a function of magnitude to take into account saturation 
effects. 

These parameters arc then used to classify all images 
into one of four categories: stellar (consistent with the 
magnitude- and position-dependent point spread function, 
cl= —1), non-stellar (a measurably extended source, cl=l), 
merged objects (sources with two local maxima within a 
single set of connected above-threshold pixels, cl=2), and 
noise (objects with nonphysical morphologies, cl=0). 

For further details of the principles involved, see Maddox 
et al. (1991a,b). Very bright images can often be misclas- 
sified, since the limited set of parameters does not provide 
an adequate description and the background-following al- 
gorithm attempts to track over them in order to detect 
the faint images in the source halos. The merged/non- 
stellar boundary is not as reliable as the stellar/non- stellar 
boundary, so merged stars are often found in the non- 
stellar list (with a smaller number of galax;ies in the merged 
list). Some objects classified as noise are real; objects 
found on both plates are the obvious examples. Objects 
classified as noise which match FIRST sources arc also 
likely to be real, and we do not generally exclude these 
from our analysis. 

Bright objects {e.g., 0,E < 13) cover a large num- 
ber of pixels in the APM scans and, as a consequence, 
magnitude and source-size estimates arc very sensitive to 
small uncertainties in the plate sky level and details of the 
background-following algorithm; as a result, large uncer- 
tainties in the parameter estimations can result, and very 
bright sources can even be completely missing from the 
catalog. In addition, bright galaxies with complex surface 
brightness distributions can be broken up into a swarm 
of discrete sources. At fainter magnitudes, the limitations 
of the plate material make reliable separation of stellar 
and non-stellar sources problematic. Since our goal here 
is to determine the completeness and reliability for the 
catalog of radio source counterparts, we do not attempt a 
comprehensive analysis of the APM catalog's classification 



accuracy. Instead we perform several straightforward com- 
parisons with existing catalogs of bright stars and bright 
galaxies and with CCD images which serve to character- 
ize the APM completeness and reliability with respect to 
FIRST identifications. 

3.6. APM Reliability at Bright Magnitudes 

3.6.1. Galaxies 

We have examined the DSS images, APM O- and E- 
platc catalogs, and FIRST images in the vicinity of 1000 
entries in the UGC galaxy compilation (Nilson 1973). The 
UGC purports to be an angular-size limited sample for ob- 
jects with blue diameters > 1'; it is complete for V < 14.5, 
but includes entries down to V = 18. This provides a di- 
rect measure of the completeness of our identification of 
FIRST sources with bright galaxies and the accuracy with 
which the galaxy parameters are reported in the APM cat- 
alog. 

Of the 1000 UGC catalog entries, seven were found to 
have no galaxy brighter than V = 17 within 3' on the 
DSS. Another dozen entries had only very small-diameter 
objects (< 20") in the vicinity, while three were extended, 
but very low surface brightness objects. Finally, in seven- 
teen cases, there were two (and in one case three) galaxies 
of comparable brightness within 1-2' of the UGC posi- 
tion, and it was not possible without further investigation 
to distinguish which was the cataloged object. 

In only 1 of the 1000 cases was no object present in 
the APM catalog, and that case corresponded to one of 
the extremely low-surface brightness galaxies; in virtually 
all cases, in fact, the UGC galaxy was detected on both 
plates. Over 275 of the 1000 UGC galaxies arc detected 
in the FIRST catalog, and all of these correspond to ob- 
jects in the APM source list. Thus, bright galaxies missing 
from the APM catalog are not a source of incompleteness 
in identifying FIRST counterparts. 

Classifications, colors, and magnitudes for bright APM 
galaxies are, however, more problematic. First, for the 
reasons discussed above, colors and magnitudes for bright 
galaxies are at best estimates, and can be grossly in er- 
ror for the brightest objects (E < 12; see §3.3). Second, 
of the 1000 UGC galaxies, 28 were classified as stellar on 
both plates; 26 were listed as stellar on the blue plate 
while correctly classified on the red, while 82 were classi- 
fied as stellar on the red plate, but correctly identified as 
extended on the blue plate. Among the i^ZRS'T-detected 
galaxies, stellar misclassifications had a higher frequency 
by a factor of 2 to 3 (~ 20% in total), presumably because 
many of these radio detections represent galaxies with a 
bright active (stellar) nucleus. 

For ~ 80 galaxies, entries appear in both the O and 
E catalogs, but their centroids are sufficiently far apart 
that they are not identified as the same object, and so re- 
main listed as separate objects (without measured colors) 
in the merged catalog. These arc slightly over-represented 
among the radio detections, most likely as a result of the 
high radio-detected fraction of interacting/merging galax- 
ies whose complicated surface brightness profiles confuse 
the source-finding algorithm. Finally, somewhat less than 
4% of the UGC galaxies have size and shape parameters 
significantly discrepant from the images, either as a conse- 
quence of including a nearby star in the profile, breaking 
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Fig. 1. — The accuracy of APM magnitudes as determined from a comparison of stellar objects taken from two different plates in the 
narrow overlap regions near the POSS-I plate edges, (a) Magnitude differences using the original APM magnitudes. The dashed line is a 
Gaussian with cr = 0.45. (b) Magnitude differences after photometric calibration using the APS catalog. The dashed line is a Gaussian with 
<T = 0.30. The recalibrated magnitudes are substantially improved. 




Fig. 2. — The differences between the calibrated and raw APM magnitudes as a function of magnitude for the (a) E and (b) O plates 
derived from the 148 POSS-I plates covering the FIRST survey area. 
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up a large galaxy into many components, or some other 
such error; these are also more common (6%) among the 
radio detections owing to the high fraction of interacting 
systems FIRST detects. 

In summary, then, the APM catalog is > 99% complete 
for bright galaxies in the sense that it contains at least one 
entry at the galaxy location, and in ~ 97% of the cases, 
the object is classified as a galaxy on at least one of the two 
plates. The descriptions of the magnitudes, shapes, and 
colors of the optical objects, however, arc; subject to large 
errors, and should be checked by examination of the DSS 
and/or other catalog resources such as NED before they 
are used for purposes other than simple identifications. 

3.6.2. Stars 

As noted in the Introduction, the number of stars de- 
tected at centimeter wavelengths is small, and the number 
of quasars brighter than = 14 is even smaller. Thus, the 
completeness of the APM for bright stellar objects is not 
a major issue in the identification of FIRST counterparts. 
Nonetheless, for the record, we briefly comment on APM 
reliability and completeness for bright stellar objects. 

Helfand et al. (1999) report on the detection of 21 ra- 
dio stars in the FIRST region under discussion here with 
magnitudes in the range 1.6 < V < 14.1. All but the 
three brightest objects (V < 4.5) arc recorded in the APM 
catalog, and more than half are classified as stellar on at 
least one plate. Since most true stellar identifications with 
radio sources are impossible without proper motion infor- 
mation, however, and the number of extragalactic bright 
stellar counterparts is vanishing small, any incompleteness 
at bright magnitudes for stellar objects in the APM cata- 
log is largely irrelevant to FIRST source identification. 

3.7. APM Reliability at Fainter Magnitudes 

In order to assess the completeness and reliability of 
the APM counterparts list in the magnitude range 15 < 
E < 20, we use statistics from work in progress to identify 
FIRST sources in the 16 deg^ deep I-band imaging sur- 
vey (/ < 24) of Postman et al. (1998b; see Helfand et al. 
1998 for a preliminary report on this FIRST identification 
program) . 

A total of 345 counterparts to FIRST sources discovered 
in the I-band CCD data were selected for comparison with 
the APM catalog; this included all 323 objects with / < 20 
and the additional 22 stellar objects with 20 < / < 21. All 
123 objects with / < 18.0 were detected on both plates 
and over 90% (75%) are correctly classified on the E (O) 
plates. Over 93% of the 87 objects with 18 < / < 19 were 
also detected on both plates, although the classifier per- 
forms somewhat less well, with 63% correctly classified on 
the E plate. For / > 19, more than 40% of the /-band ob- 
jects are still detected on the E and/or O plates, of which 
half are correctly classified on the E plate. This near to 
the plate limit, it is clearly difficult to distinguish between 
stellar and extended objects; nonetheless, only one of the 
27 spectroscopically identified quasars with 17 < 7 < 21 
in this /-band field is misclassified as non-stellar on both 
plates. Similarly, we have found that fewer than 5% of the 
non-Seyfert quasi-stellar objects in the Veron-Cetty and 
Veron (1998) catalog with R < 18 are mistakenly classi- 
fied as galaxies on both plates (Gregg et al. 1996). 



Using an APM magnitude-limited sample at E < 19.0, 
we find that 62% of galaxies are classified correctly on both 
plates. An additional 28% are classified correctly on the 
E plate and are at or below the plate limit on the O plate, 
while another 4% arc classified correctly in O and incor- 
rectly in E; only 6.5% are misclassified on both plates. For 
stellar images, the score is similar, with 14/26 correct on 
both plates and another seven classified correctly on one 
plate. In summary, for E < 19, 92% of all objects are 
classified correctly on at least one plate. 

As a further check on the incompleteness of the APM 
catalog, we have compared a subset of the database with 
the APS catalog of the POSS-I (Pennington et al. 1993). 
The APS catalog contains only objects detected on both 
plates, so it does not include the faintest APM objects 
which are often detected on only a single plate. This is 
nonetheless a very useful test of completeness and accu- 
racy for the APM catalog, since the APS started with 
the same plate material but used completely independent 
scanning hardware and processing software. 

The APS catalog test set contains 61,000 unique objects 
that match 71,000 different APM objects within a match- 
ing radius of 10"; the objects were drawn from areas on the 
plates within 20" of FIRST souices, and thus include both 
real radio source identifications and random background 
sources. There are 532 APS objects (0.9%) that have no 
matching APM source. Most of these are very near the 
catalog limit: 364 have E > 20 or O > 21, leaving only 
168 (0.3%) that are reasonably bright. Checking the Dig- 
itized Sky Survey reveals that the great majority of these 
objects are blended with other nearby objects in the APM 
catalog (and usually classified as such), leading to poor 
agreement in the positions from the two catalogs. Fur- 
thermore, we should note that we have not examined the 
plates by eye at the locations of these "missing" sources, 
and some could be spurious APS entries which arc legit- 
imately absent in the APM catalog. Thus, the fraction 
of objects that is simply missing from the APM catalog 
above i5 = 20 (O = 21) is very small, certainly less than 
0.1%. 

Occasionally the positions measured by the APM on the 
O and E plates are sufficiently different that the red and 
blue detections are not recognized as being a single source 
and so two entries (one E-only and one O-only) appear 
in the APM catalog (§3.6.1). This can lead to incom- 
plete radio source identification. In our APS test sample, 
there are ~ 730 cases (1.2%) where a single APS source 
matches a close pair (< 10") of £^-only and O-only ob- 
jects from the same APM plate. We consider these objects 
to be cases where a single source has been split into two 
catalog entries. The majority of these split objects are 
bright (Fig. 3) and are further examples of the complex 
and/or blended objects that occasionally caused trouble 
for the bright UGC galaxies discussed above. The median 
E magnitude of these objects is 16.1; the fraction of such 
objects is < 1% for the magnitude range 15-20 in which 
most radio counterparts are found. 

4. FIRST AS AN ASTROMETRIC STANDARD 

4.1. Introduction 

As we show in WBHG, all unresolved FIRST sources 
down to the survey limit have 90% confidence positional 
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Fig. 3. — The distribution as a function of APM E magnitude of the fraction of objects that appear in the APM catalog as a pair of sources, 
one detected only on the O plate and one detected only on the E plate. These split objects are quite bright and arc typically cither large 
complex galaxies or close, blended objects. The fraction of such objects in the range of most FIRST identifications (15 < E < 20) is very low. 



uncertainties of < 1"; furthermore, the absolute reference 
frame is tied to the VLB reference frame with an offset of 
< 0.03" and an rms uncertainty of < 0.15" in each coor- 
dinate. This level of precision is unprecedented for such a 
high-surface-density catalog {e.g., the Guide Star Catalog 
entries have an rms uncertainty of 0.4" with systematic 
errors of up to 1"), making FIRST a useful astrometric 
calibration standard for any astronomical catalog which 
contains a significant number of FIRST source counter- 
parts. In this section, we use the radio source positions to 
calibrate and correct the APM intraplate astrometry er- 
rors and the plate-to-plate shifts, allowing an improvement 
by nearly an order of magnitude in POSS-I positions. 

The FIRST catalog contains about 90 sources deg~^ 
of which ~15% (see §5.2) have APM optical counterparts 
within 2", yielding -^600 matches per POSS-I plate. This, 
in effect, provides a dense grid of astrometric standards 
across each plate. More important, the optical counter- 
parts to FIRST somces are generally faint, and thus com- 
plement the astrometric calibration that can be achieved 
using bright standards such as the PPM and Tycho cat- 
alogs that have been used heretofore as the basis of the 
APM plate solutions. 

4.2. Overall APM Astrometric Offset 

Whilst working on radio identifications for the Jodrell 
Bank-VLA Astrometric Survey (JVAS - Patnaik et al. 
1992; Hook et al. 1996; Snellen et al. 2001), it was noted 
that there was evidence for a shift between the mean Right 
Ascension and Declination of the APM reference frame 
and the VLA reference frame. The mean shift was derived 
by taking the median of the nearest optical counterparts 
to the survey's flat-spectrum radio emitters and yielded: 

• RA^PM (corrected) = RA^pm (original) + 0.35 
±0.05 arcsec 



• DecAPM (corrected) = Dccapm (original) -|- 0.35 
±0.05 arcsec 

We verified that this effect was not intrinsic to the PPM 
reference frame used for the APM plate solutions at that 
time by comparing the optical positions of a set of VLBI 
Radio Reference Frame objects (Johnston et al. 1995) with 
optical counterparts on APM scans and found that the 
same systematic shift was present. We then remeasured 
POSS field 1393 centered on the North Galactic Pole in 
two orientations: with the normal scanning direction and 
with the plate rotated by 180 degrees. We found that, after 
application of the original astrometric alignment using the 
PPM stars, the derived positions for faint (> 14*'' magni- 
tude) objects differed by ~ 0.7" in both coordinates. This 
indicated that the APM measurement system was intro- 
ducing a systematic shift in the positions of bright objects 
with respect to fainter ones. The origin of this effect is still 
under investigation. Our first step in comparing FIRST 
sources and the APM catalog, then, was to take out this 
systematic shift of 0.35" in both coordinates. 

4.3. Intraplate Errors 

The original POSS-I plates suffer from various distor- 
tions resulting from the stress induced by the plate holder, 
the vignetting of the 48-inch Schmidt telescope, etc., which 
are reproducible from plate to plate (Irwin 1994). In Fig- 
ure 4(a), we display the pattern of residuals from the 
TYCHO reference stars derived by stacking the APM- 
TYCHO offsets for all such stars on 148 plates. Systematic 
errors reach 1.3" in some regions of the plate. A correc- 
tion for this effect has, as noted above, been applied to 
all APM catalog positions using a correction map derived 
from the Tycho catalog. 

We compared these corrected APM positions with those 
of FIRST radio sources over the 148 POSS-I plates which 
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Fig. 4. — (a) The distortion in raw APM catalog positions as a function of position on POSS-I plates derived by comparison to reference 
positions from the TYCHO catalog of astromctric reference stars. The lines show the average position shifts in 20' X 20' regions across the 
plates (6667 pixels = 1 degree); the scale of the shifts is shown on the right. There are large, systematic distortions that are not modeled 
by the APM plate solution, (b) The distortion remaining in the APM positions after correction for the distortions shown in (a). These 
are derived by comparing corrected APM positions to FIRST positions for APM-FIRST matches. The distortions are small near the plate 
center, but beyond ^ 2.5° from plate center there is a strong, increasing radial distortion. This is a magnitude-dependent position error in 
the APM catalog; the positions of bright sources (such as TYCHO stars) are properly corrected by application of the distortion map in (a), 
but positions of faint objects (such as the optical counterparts to FIRST radio sources) are actually made worse by the bright-star correction 
map. 
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Fig. 5. — The radial oflFset between APM and FIRST positions as a function of distance from plate center; each point represents one of the 
20' X 20' grid cells into which we have divided the plates. 
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subsume the current survey region. Including all sources 
with optical counterparts having offsets < 2.0" (~60,000 
objects with a false rate of ~ 5% — see §5) yields the re- 
sult displayed in Figure 4(b). Residuals within a radius of 
2.5 deg of the plate centers are substantially reduced com- 
pared with Fig. 4(a), but large radial errors with mag- 
nitudes up to 1.2" are now seen around the plate edges 
(Fig. 5). We believe this effect is a consequence of a com- 
bination of telescope vignetting and plate saturation of 
the bright stars in the Tycho grid: vignetting distorts the 
stellar images in the radial direction, producing elliptical 
isophotes whose centroids are systematically shifted to- 
ward the plate center. The FIRST counterparts, largely 
much fainter objects for which saturation is not a problem, 
reveal this distortion. 

To correct for this effect, we have chosen to calibrate 
the errors using all FIRST counterparts with 15.0 < E < 
19.5; the bright limit is chosen such that the magnitude- 
dependent errors described above will be unimportant 
(and only excludes 6% of the matches), while the faint 
limit is imposed in order to include only those objects 
for which the random errors are small. The correction 
proceeds iteratively. We calculate a mean offset (in both 
the radial and tangential directions) for all sources in the 
stacked image in cells 20' on a side (20x20 cells cover 
the plate). We then apply the derived offsets to all APM 
objects in the catalog based on their x-y plate positions, 
rematch to FIRST, and then recalculate the mean offset 
for the cell. The process is repeated until it converges 
(~4 steps). An iterative approach is required because 
the APM position shifts are comparable to the 2" APM- 
FIRST matching radius. The resulting APM-FIRST po- 
sition errors are much improved: the maximum residual 
mean offsets are less than 0.18" with an rms of only 0.06" 
over the whole plate. 

This intraplate correction is strictly valid only for the 
range of magnitudes used in the calibration (15.0 < E < 
19.5). An examination of the offsets for fainter coun- 
terparts suggests that the corrections work well to the 
plate limit, although the positional uncertainties increase 
slightly owing to the statistical fluctuations which are in- 
evitable in deriving positions for objects near the detection 
threshold. For the brighter objects, however, the derived 
positions of candidate counterparts to FIRST sources suf- 
fer from the same type of distortion evident in the Tycho- 
based solution, presumably in a magnitude-dependent 
way. Rather than attempting to derive the functional 
form of the magnitude dependence from the small num- 
ber of bright i^/i?5'T counterparts, we recommend that for 
very bright objects either the separation for an acceptable 
match be increased or the offset be computed using the 
APM astrometry without the FIRST astrometric correc- 
tion. For example, in the FIRST Bright Quasar Survey 
(White et al. 2000) we include all sources with optical 
matches closer than 1.2" in eit/ier astrometric frame, a pro- 
cedure which adds 24 objects to the 1214 selected using the 
FIRST astrometric solution alone. Since the typical dis- 
crepancy between the two solutions is < 1" and our mini- 
mum matching criterion is > 1", no genuine matches will 
be missed if this procedure is followed. However, since the 

^ We display only the region within the 1150x1550 pixel image from 
will be extracted from the adjacent coadded images. 



overwhelming majority of optical counterparts are faint, 
we have used only the matches with the FIRST correc- 
tions applied when compiling the match statistics in this 
paper. 

This use of FIRST sources as astrometric calibrators 
assumes, of course, that there are no systematic errors in 
the radio positions. In Figure 6, we display the analo- 
gous plot to Figure 4(b) for the FIRST catalog: the x-y 
positions of each FIRST source in the radio "plate" {i.e., 
coadded image — see BWH) arc extracted and the; mc;an 
offset from the APM matches to these sources (binned in 
2' cells) is calculated. The uniformity is excellent and the 
magnitude of the errors is very small: they are all < 0.13" 
with an rms of 0.05" over the whole field'^. This result, 
coupled with our global astrometric calibration described 
in WBHG suggests our confidence in FIRST astrometry 
is justified, and that our radio catalog can be utilized rou- 
tinely to aid the astrometric calibration of other catalogs. 

4.4. Plate-to-plate Shifts 

Having established the global APM offset and calibrated 
out the reproducible intraplate errors, we are left only with 
translations of individual plates with respect to our radio 
reference frame. The mean FIRST-APM offsets for each of 
the 148 plates in the survey region are shown in Figure 7. 
These offsets were computed using an iterative approach 
similar to that used to derive the plate distortion. There 
is no obvious trend with Right Ascension or Declination, 
but shifts remain at levels ranging from 0.01" to 0.70"; the 
distribution of errors is roughly Gaussian in both coordi- 
nates with an rms of ~ 0.16" and an apparent tendency for 
large offsets in one coordinate to be matched by those in 
the other. We remove these shifts by simply subtracting 
the shift computed for a plate from the APM positions 
of all objects on the plate. Note that for plates which 
fall completely within the boundaries of the FIRST sur- 
vey, this calibration is final; however, for those plates only 
partially covered by the existing survey, the global offset 
will be redetermined as new data accumulate. The size 
of any subsequent correction is expected to be no more 
than ^0.1". The current values, along with the intraplate 
correction matrix are available on the FIRST Website. 

After correcting for the plate shifts, we have recom- 
puted the plate-based distortion maps (Fig. 4) to deter- 
mine whether the systematic plate shifts affect them. The 
changes in the distortions are small (< 0.1"). 

Having completed the astrometric calibration, we can 
now proceed to generate a catalog of optical identifications 
for the FIRST survey. 

5. OPTICAL IDENTIFICATION OF FIRST RADIO SOURCES 

The primary condition for assigning an optical counter- 
part to a radio source has traditionally been positional co- 
incidence; while other factors can enhance or detract from 
the probability that the association is correct, these are 
generally second-order effects. If the positional accuracies 
for both the radio and optical catalogs are very high, posi- 
tional coincidence alone can afford a high degree of confi- 
dence. In BWH, we calculated the expected error rate for 
associations between the APM and several hypothetical 

h sources actually enter the catalog; sources closer to the boundaries 
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Fig. 6. — The distortion in FIRST catalog positions as a function of position within the FIRST coadded images. The scale of the vectors 
is shown on the right; 1 pixel = 1.8". No systematic errors in the FIRST positions are evident, indicating that we have successfully corrected 
the geometric distortions in the FIRST images. 




Fig. 7. — The mean shifts for each APM plate in the FIRST survey area derived by comparing APM positions (corrected for the intraplate 
distortions in Fig. 4) to FIRST positions. Each point shows the RA and Dec shifts for a plate with Icr error bars. 
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VLA radio surveys with differing angular resolutions. The 
reliability of associations based on positional coincidence 
were presented as a function of optical source density on 
the sky, and as a function of apparent V magnitude at 
the North Galactic Pole. The results demonstrated the 
strong dependence of both reliability and completeness on 
the angular resolution of the radio survey. 

In reality, the actual error rate in radio/optical associ- 
ations is a complex problem depending on the details of 
the relative astromctric accuracy of the two data sets, the 
optical morphology of candidate objects, and radio source 
morphology. In this section we quantify empirically the re- 
liability of associations between FIRST radio sources and 
optical counterparts on the POSS-I plates as measured by 
the APM. 

5.1. Background Rates Using Closest Matches 

In our analysis of the FIRST-APM match rate, we re- 
tain only the closest optical match to each radio source 
and vice versa. While the optical catalog we present below 
includes all matches out to 20", not just the closest one, 
using only the closest matches when describing the iden- 
tification statistics makes understanding the background 
and true match rates easier. We want to know the num- 
ber of radio sources that have optical counterparts, but 
we want to avoid double-counting {e.g., a particular opti- 
cal object may be the closest coimterpart for two differ- 
ent radio sources.) Also, neither catalog can have pairs 
of sources with separations of less than a few arcseconds, 
which means that there is effectively a "hole" in the cat- 
alog around each radio or optical source. This makes it 
very difficult to calculate chance coincidence rates, since 
true matches suppress additional false matches in a com- 
plicated manner. 

We can avoid both of these problems by using only the 
closest matches. This also modifies the chance background 
rate, since once a source has a counterpart there can be no 
additional counterparts at wider separations. This effect 
is straightforward to understand and incorporate into the 
background computation: the area being searched for new 
matches at a given radius is contributed only by FIRST 
sources that do not already have a closer optical counter- 
part. Another complication when using closest matches is 
that the matching circles for close pairs of radio sources 
can overlap. An optical source that falls in the overlap re- 
gion between two radio sources can only be assigned as the 
counterpart to one of them. This is also a mathematically 
tractable problem. Both effects are treated as a reduction 
in the effective sky area that is searched as a function of 
matching separation. 

When both of these effects are taken into account, it 
is possible to predict the number of chance (false) coin- 
cidences as a function of the radio-optical separation. A 
simple model for this prediction that assumes a uniform, 
random distribution of optical sources on the sky predicts 
more coincidences at large separations than are actually 
seen, owing to the variation in optical source density Popt 
on the sky. In the presence of a varying optical source 
density, radio sources that remain unmatched at large sep- 
arations are more likely to be found in low-density areas 
of the sky; consequently, the effective mean optical source 
density decreases as the separation increases. If the vari- 



ation 5 Popt/ Popt is not too large, this can also be modeled 
fairly simply. For an ensemble of N radio sources, the 
effective background source density is 

MAr) = -popt-alJ^^^^^) , (6) 



where popt is the mean source density, Uopi is the rms vari- 
ation in the density, Ai{r) is the area already searched for 
source i, and (pi{r) is the angular part of the matching 
circle being searched for source i. For an isolated FIRST 
source with no nearby neighbors (so there are no overlap- 
ping error circles) and an optical match at radius ro, we 
have 

?!>i = 27r , r<ro (7) 
= , r>ro , (8) 

and 

= Trr^ , r< ro (9) 

= 7rro , r > ro . (10) 

For "gregarious" FIRST sources in close groups, the 
matching circles overlap. Then <pi will include only the 
portion of the error circle at r that would not be assigned 
to another FIRST source, so that in general (pi < 27r. 

Our background model, then, has two parameters: the 
mean source density popt and the variance in the density 
"'opt- Values for these quantities were estimated from a 
spurious match catalog generated by offsetting the coordi- 
nates of the radio sources by 5' to the south. We use this 
procedure instead of simply adopting the match rate be- 
tween, say, 10" and 20" around the real source positions in 
order to minimize any enhancement in the false rate result- 
ing from the presence of real radio-optical associations at 
large separations. Such associations result both from op- 
tical counterparts to multiple-component FIRST objects 
(where the optical position need not match any one of 
the cataloged FIRST components closely) and clustering 
of galaxies around FIRST objects (which will produce an 
uncharacteristically high optical source surface densities in 
the vicinity of the radio sources.) The parameters for the 
background due to various types of APM objects are given 
in Table 1. 

5.2. Match Results 

In Figure 8, we display the result of matching all 382,892 
radio sources in our catalog to the astromctrically cor- 
rected APM catalog. The plot shows the cumulative ex- 
cess of matched sources over the background of chance 
coincidences as a function of the offset between the ra- 
dio and optical positions. We estimate that 98% of the 
APM sources within 1" of a FIRST source are physically 
associated with the radio source; 42,400 sources meet this 
criterion. Even out to 2", 94.5% of the 59,700 associations 
are real. Figure 8 indicates that some real matches occur 
out to > 10" although the reliability decreases steadily as 
the separation increases. Integrating under the curve of 
Figure 8 out to 4" implies 61,800 real associations (16% of 
all radio catalog entries). Inside the 1" radius, 16% of the 
optical counterparts are classified as stellar on both plates, 
while 41% are classified as non-stellar on both plates; in 
the remaining cases, the classifiers disagree or the object 
is only detected on one plate. Note that, at faint magni- 
tudes, classification becomes difficult owing to the limited 
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number of pixels above threshold; in addition, active nu- 
clei can make a galaxy appear stellar in the blue band, and 
the bulges of faint ellipticals and SOs are generally unre- 
solved. Thus, galaxies are increasingly classified as stellar 
as the plate limit is approached. 

In Table 3 we tabulate the completeness and reliability 
of objects that are only detected on a single red or blue 
plate. These have been estimated by shifting the radio po- 
sitions by 3 arcmin in declination and running the match- 
ing analysis in an identical manner on this shifted dataset. 
This shows that 92.7% of the blue-only matches within 1" 
are real, and 97.5% of the red-only matches within 1" are 
real matches. This is a lower limit on whether these single 
band detections are real objects since we expect some of 
the the chance associations to be real celestial objects. 

The matching results are more complicated to analyze 
when we take into account the complex radio sources that 
get broken into two or more components in the FIRST cat- 
alog. We call these "gregarious" sources (as distinguished 
from isolated sources) because they are found in clumps on 
the sky. Note that some gregarious sources are simply the 
chance superposition of two unrelated radio emitters at 
different distances; nonetheless, from our vantage point, 
they are gregarious, and confuse the matching statistics 
in a manner similar to that of the real multi-component 
objects. We define the sociological boundary between gre- 
garious and isolated at 60"; i.e., any source with no other 
catalog entry within a radius of 60" is classified as isolated, 
and all other sources are labeled gregarious. Such a fixed 
boundary is arbitrary, and, indeed, a small number of very 
extended objects will be incorrectly classified as isolated; 
we discuss this matter further in §5. 

In Figure 9, we display the number of isolated matched 
sources as a function of the offset between the radio and 
optical positions in 0.1" bins, normalized by the annular 
area of each bin. Figure 10 displays the ratio of the mrni- 
ber of matches to the predicted false rate as a function of 
separation. 

Both the background rate and the typical angular 
separation between FIRST and APM positions depend 
strongly on the optical classification. Figure 11 shows the 
distribution of separations for sources classified as stellar 
on both plates or non-stellar on both plates. The FIRST- 
APM positions typically differ by more for galaxies, which 
have less well-determined optical positions. However, it 
is easier to identify galaxies confidently as FIRST coun- 
terparts because the background rate for galaxies is four 
times smaller than for stars (see Table 1). 

We quantify this effect by calculating the angular radius 
which contains 90% of the real associations as well as the 
percentage of false matches within the same radius. For 
example, these numbers are 4.2" and 17% for the distri- 
bution of all FIRST/ APM sources as shown in Figure 8. 
If we restrict the match to APM objects that are classified 
as galaxies on both plates (Fig. lib), the values are 4.0" 
and 2.7%. The 90% radius is quite similar to that for all 
associations, since galaxies constitute the majority of all 
identifications; however, the reliability for galaxies is much 
higher because the background rate for galaxies is smaller 
(Fig. 11). By comparison, if we restrict ourselves to APM 
sources that are stellar on both POSS plates (Fig. 11a), the 
90% radius and error rate are 1.7" and 6%. The smaller 



footprint of stars on the POSS-I leads to more accurate 
optical positions (hence the smaller 90% radius). These 
and other cases are summarized in Table 2 and in Fig- 
ure 12, which displays the completeness and reliability of 
matches as a function of separation. 

5.3. The Effects of Radio Morphology 

The statistical agreement between FIRST and APM 
positions is also affected by the morphology of FIRST 
sources. If we restrict the FIRST sample to point-like 
radio sources (but include all APM sources), the 90% as- 
sociation radius is 2.0" with a false rate of 5%. Limiting 
the discussion to sources that are point-like in both the 
radio and the optical results in 90% of the matches within 
1.1" with a 2.4% false rate. 

Figure 13 compares the histogram of separations for 
isolated sources with that for the gregarious sources. It 
can be seen that gregarious sources contribute most of the 
matches at large separations; indeed, they show evidence 
for a statistically significant excess of matches even beyond 
15". Many of the gregarious FIRST sources are compo- 
nents of classical double radio sources, which may have no 
radio component at the center of the double. Consequently 
the optical counterpart will typically be ^ 1/2 the double 
separation from each component. We now begin to explore 
counterpart identification strategies for such sources. 

5.4. Optical Counterparts to Double Radio Sources 

The wide variety of radio source morphologies (not to 
mention the extent to which their appearance depends on 
the resolution of the survey) makes it difficult to design a 
robust algorithm that predicts the locations of the optical 
counterparts to multiple-component sources. We present 
here an empirical approach to the problem for the simplest 
class of such sources - isolated doubles. 

Figure 14 displays the distribution of ~ 45, 000 isolated 
pairs of FIRST sources as a function of their separation. 
The minimum value of 3" is imposed by the source detec- 
tion algorithm (WBGH). The fall-off in the number of de- 
tected doubles beyond ^ 10" reflects a number of factors: 

1) the 5" survey beam resolves out very extended features 
and will therefore not detect sources on the largest scales, 

2) the exclusion of triples and other multiple component 
sources from this subsample will preferentially remove ex- 
tended sources on larger scales, and 3) there is a real de- 
cline in the number of large angular-diameter objects. For 
comparison with the APM catalog, we have selected ob- 
jects with separations in the range 8" < d < 30". 

The lower portion of Figure 15 displays the distribution 
of optical objects in the vicinity of the 21,579 pairs of radio 
sources meeting this criterion. The x-axis for each pair is 
defined by the line joining the centroids of the two objects; 
the scale is normalized to the separation, with the origin 
chosen as the brighter of the two components. The upper 
panel shows a histogram for all optical objects within 1.5" 
of the line joining to two components; the expected false 
rate from a uniform distribution of sources with the same 
mean surface density as the FIRST snrycy is shown as the 
dashed line. Several features of the distribution are im- 
mediately apparent. There is a large concentration of op- 
tical objects coincident with the brighter of the two radio 
components, suggesting a core-jet morphology. A roughly 



14 



McMahon, White, Becker, & Helfand 



8-10 



m 6-10^- 



E 4-10^- 



E 

o 



2-10^ - 




5 10 
Separation (arcsec) 



Fig. 8. — The cumulative number of FIRST-APM matches as a function of separation. Only the closest APM object to each FIB.ST 
source is included; both isolated and "gregarious" FIRST sources (those with near neighbors) arc included. The background coincidence rate 
(dotted line) has been subtracted. The inset shows the cumulative fraction of sources as a function of radius with an expanded scale. The 
vast majority of matches within 2" are clearly true identifications; a few true identifications are found out to > 10" separations. A total of 
~ 73, 000 FIRST sources have optical counterparts in the APM catalog. 
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Fig. 9. — A histogram of the number of matches as a function of separation between the FIRST and APM positions. Only the closest APM 
match to each FIRST source is included, and only isolated FIRST sources have been included to avoid the ambiguities in the background 
computation for multiple-component FIRST matches (see Fig. 13.) The histogram has been normalized by the product of the area of each 
annular bin and the number of FIRST sources, so it gives the mean number of APM matches per unit area for each FIRST source. The 
background rate of coincidental matches is shown in more detail in the inset (error bars are Icr and the dashed line shows the expected chance 
coincidence rate.) 
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Fig. 10. — The ratio of the number of APM-F/i?5r matches to the expected number of coincidental matches as a function of separation. As 
for Fig. 9, only the closest APM matches to isolated FIRST sources arc included. This function shows the effect of increasing the matching 
radius between the catalogs. At a separation of ~ 2.5", the background rate per unit area is approximately equal to the true match rate, so a 
small increase in the matching circle size leads to about equal numbers of true and false matches being included. Within 2" the vast majority 
of the matches (98%) are real associations; when non-isolated FIRST sources are included, this fraction declines to 95%. 
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Fig. 11. — The match rate as a function of separation for (a) APM objects classified as stellar on both the POSS-I plates and (b) APM 
objects classified as iioii-stcUar on both plates. As for Fig. 9, only the closest APM matches to isolated FIRST sources are included, and the 
insets show the counts to larger radii with a much expanded scale. Galaxies are less concentrated to small separations than stars (mainly 
because the optical positions are less well-determined), but the background rate for galaxies is four times smaller than that for stars, which 
makes matches to galaxies reliable to larger radii than those for stars (see Fig. 12). 
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Fig. 12. — Completeness and reliability as a function of separation for (a) objects classified as stellar on both APM plates and (b) objects 
classified as non-stellar on both plates. The completeness rises almost monotonically from zero to one as the separation increases (the small 
deviations from moiiotoiiieity are due to the background subtraction). The differential reliability, which measures the probability that sources 
at a given separation are true associations rather than chance coincidences, declines toward larger separations. The solid lines show the rates 
for isolated FIRST sources, while the dashed lines are for gregarious FIRST sources. Stellar matches to isolated sources axe very closely 
concentrated toward small separations, and few stellar matches beyond a few axcseconds are real. On the other hand, because the background 
rate for non-stellar matches is much smaller (see Fig. 11), matches with galaxies are reliable to much larger radii. 
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Fig. 13. — The normalized histogram of separations for isolated FIRST sources (with no other cataloged FIRST objects within 60") and 
"gregarious" FIRST sources (which do have neighbors in the catalog.) The background coincidence rate has been subtracted. Clearly most of 
the matches at large separations are attributable to the multiple-component gregarious sources, where the position of the optical counterpart 
is not as tightly coupled to the FIRST radio position as it is for isolated sources. Developing a matching strategy for gregarious FIRST 
sources and determining the chance coincidence rate is a difficult problem which we begin addressing in §5.4. 
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Fig. 14. — The distribution of angular separations of all isolated pairs of FIRST sources (no third component within 120"). The minimum 
separation of 3" is set by the resolution of the FIRST images and our detection algorithm; the fall-off beyond 10" arises from a variety of 
effects (see text). The dashed line represents an estimate of the chance coincidence rate. At a separation of 40", approximately half the 
doubles are physically associated and half are chance superpositions of unrelated sources 



equal number of identifications is found approximately half 
way between the two components with a slight bias in the 
direction of the brighter lobe; the broader spread in the y- 
direction for these counterparts reflects in part the bending 
of radio lobes as a consequence of their interaction with 
the intergalactic medium {e.g., Blanton et al. 2000, 2001). 
Finally, a much smaller fraction of the identifications is 
coincident with the weaker of the two lobes. The over- 
all identification rate derived from integrating the upper 
curve and subtracting the background rate is ~ 19%, sim- 
ilar to that for the radio sample as a whole. 

Dividing the optical counterparts by morphological class 
yields additional interesting results. The objects classi- 
fied as galaxies on both plates (Figure 16a) produce the 
most counterparts midway between the components, while 
the stellar objects are predominantly associated with the 
brighter radio component (Figure 16b). Figure 17 shows 
that many of these core components are small in size (panel 
a) and are much brighter (panel b) than the other compo- 
nent: the median flux ratio for sources where the optical 
counterpart matches the brighter radio component is 3.0, 
while the median for the central-component matches is 1.4. 
When the brighter of the radio components identifies the 
optical counterpart, it also tends to be the smaller of the 
two components (panel c). Likewise, for cases in which the 
dimmer source has the optical counterpart, it tends to be 
smaller. In contrast, the sources with identifications that 
lie between the two radio lobes have very narrow ranges 
of flux density and size ratios near unity, and are virtually 
all resolved. 

Using these trends in radio source component flux den- 
sity and size ratios, total extents, and optical morphology 
(in addition to magnitude and color, perhaps) it would be 

Note this number is different from that in Figure 15 because the latter 
the optical coincidence rates. 



possible to develop a reasonably reliable identification al- 
gorithm for double and multiple FIRST sources, although 
we regard such an efi'ort as beyond the scope of this pa- 
per. What is clear from the foregoing analysis is that, 
since the identification rate for complex sources is similar 
to that for single-component objects, the effect on the final 
fraction of radio emitters identifiable at the POSS-I plate 
limit simply scales with the number of radio components 
removed from consideration as a consequence of their as- 
sociation with another catalog entry. Using an algorithm 
that assigns a probability to component associations based 
on flux density and separation, we find that, for separa- 
tions up to 120" there are 32,312 doubles^, 12,802 triples, 
and 5,716 groups of four or more sources with a probabil- 
ity of real association > 90%. This reduces our catalog 
to ~ 306, 000 discrete radio emitters and yields an overall 
identified fraction at the POSS-I limit of - 24%. 

6. CHARACTERISTICS OF THE OPTICAL COUNTERPARTS 
TO FIRST SOURCES 

The objects in the FIRST catalog span over four orders 
of magnitude in radio flux density, and the identiflcation 
rate as a function of flux density is displayed in Figure 18. 
Even for sources as bright as 1 Jy, only one-third of the ob- 
jects have counterparts brighter than the plate limit. The 
identification rate falls monotonically with flux density to 
20 mJy where it reaches a minimum of ~ 12% and then 
climbs again toward the survey limit. The rise at low flux 
densities reflects the change in the composition of the ra- 
dio source population which is dominated by distant AGN 
at high flux densities and by nearby star-forming galaxies 
below 1 mJy (Condon 1992). The fraction of sources classi- 
fled as stellar on at least one of the two POSS plates is also 

includes only isolated doubles so as to examine as cleanly as possible 
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Fig. 15. — The distribution of the 8220 optical objects in the vicinity of the 21,579 isolated doubles with separations between 8" and 30". 
The axes are defined such that the line joining the two component centroids is the x-axis and the normal to this line is the y-axis. Distances 
along the x-axis are normalized to the component separation with the origin defined to coincide with the brighter of the two components; the 
y-axis is in arcscconds. The lower panel shows the distribution of all optical objects, while the upper panel displays a histogram of the 5348 
objects falling within \y\ < 1.5". The false rate is shown as the horizontal dashed line. 




Fig. 16. — Using the format of Fig. 15, we display the results for optical counterparts to radio doubles divided between (a) the 1492 objects 
classified as galaxies on both plates and (b) those 1278 objects classified as stellar on both plates. 
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Fig. 17. — (a) The angular size of the brighter radio component for the 5348 doubles that have optical counterparts within 1.5" of the 
line joining them, versus the distance along that line (normalized to the double separation as in Fig. 15). Contour linos illustrate the density 
of points in various parts of the diagram at levels of 0.01, 0.02, 0.05, 0.1, 0.2, 0.3, . . . , 0.9 of the peak. A source is clearly resolved in the 
FIRST survey if it has an angular size greater than 2.5-3". (A total of 125 objects with deconvolved sizes of 0.0" are omitted from the plot; 
54 fall near dx = 0, with 51 near dx = 1 and only 12 lying between the components.) The brighter of the two components (which defines the 
x-origin) is more often point-like, (b) The ratio of the flux densities of the radio doubles with counterparts as a function of distance along 
the line joining the two components. In cases where the optical counterparts coincide with one component or the other, the components tend 
to have very unequal flux densities, whereas for cases when the optical counterpart is between the two components, the doubles are of nearly 
equal flux density. The thick line indicates the flux-weighted centroid of the source. 
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Fig. 17. — (c) The ratio of the fitted major axes of the doubles with counterparts as a function of distance along the line joining the 
components, defined as the size of the brighter component over the size of the dimmer one. When the optical counterpart coincides with one 
component or the other, the size of that component is smaller, whereas when the counterpart is found between the two components, they 
tend to be of equal size. 
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reflective of this trend (Figure 19). Above 100 mJy over 
80% of all counterparts are stellar (quasars) , whereas this 
fraction decreases monotonically until at 1 mJy, 60% of 
the counterparts are classified as galaxies; given the strong 
tendency toward the misclassification of faint galaxies as 
stellar, the true fraction of galaxy counterparts is consid- 
erably higher. 

One way to explore the types of objects associated with 
FIRST sources is to place the optical counterparts onto a 
color - magnitude diagram. In Figures 20 and 21 we show 
the color-magnitude diagrams for the 'stars' and galaxies 
(based on a consistent APM classification from both POSS 
plates) within 1.5" and 4" match radii, respectively. For 
comparison, these figures also display similar diagrams for 
an equal number of stellar and nonstellar sources taken 
from random positions on the same plates. 

The radio counterparts classified as stellar are clearly 
much fainter and much bluer than most 'stars'; this is il- 
lustrated in Figure 22, where we plot the fraction of the 
radio-detected 'stars' as a function of location in the color- 
magnitude plane. The detected fraction peaks at ~ 10% 
for colors around O - = for 18 < < 20. These 
objects are nearly all quasars. The increasing fraction of 
detected objects with all colors fainter than 18*** magni- 
tude is largely an artifact produced by the failure of the 
classifier to distinguish stars from galaxies near the plate 
limit; in fact, most of the faint, red stellar objects are actu- 
ally galaxies, although a handful of very red objects have 
been shown to be high redshift quasars (Hook et al. 1998). 

For the nonstellar objects, it is clear that the radio de- 
tections arc biased toward brighter magnitudes, but are 
found for galaxies of all colors. In Figure 23, we show the 
detected fraction of galaxies in the color-magnitude plane. 
The roughly uniform detection rate as a function of color 
and the bias towards brighter galaxies are apparent; up 
to 10% of the brightest galaxies are detected and even for 
E>16, >2% of all galaxies have FIRST counterparts. 

7. THE CATALOG OF FIRST OPTICAL COUNTERPARTS 

Table 4 displays a sample page of the catalog of APM- 
FIRST matches which is available from the FIRST home 
page (http://sundog.stsci.edu). Columns 1 and 2 contain 
the radio source coordinates (epoch J2000.0) and are fol- 
lowed by the offsets in Right Ascension, Declination, and 
radius (all in arcsec) of the nearest optical object falling 
within 5" of the radio position. The global APM/ FIRST 
shift, the plate-to-plate shifts, and the intraplate shifts as 
determined by FIRST (§4) have all been applied to the 
raw APM catalog positions before calculating these off- 
sets. Columns 6-8 contain the optical classification, the 
psf parameter (see §3), and the magnitude of the coun- 
terpart on the red (E) plate. The optical classification 
codes are: stellar (consistent with the magnitude- and 
position-dependent point spread function, cl= —1), non- 
stellar (a measurably extended source, cl=l), merged ob- 
jects (sources with two local maxima within a single set 
of connected above-threshold pixels, cl=2), and noise (ob- 
jects with nonphysical morphologies, cl=0). For further 
details of the principles involved, see §3 and Maddox et 
al. (1991a,b). A negative magnitude indicates that the 
object was not detected in this waveband and the abso- 
lute value is a 2a upper limit. Columns 9-11 list the 



same information for the blue (O) plate. Column 12 is 
the (O — E) color derived from columns 8 and 11 unless 
the object is only detected on one plate, whereupon the 
color is designated -9.99. All raw APM magnitudes have 
been corrected for plate to plate systematic calibration er- 
rors as described in §3.4; the uncertainties (rms) in the 
magnitudes are ~ 0.2. The POSS plate number in column 
13 is followed by selected radio parameters for the FIRST 
source, including the peak and integrated flux densities, 
and the deconvolved source major axis, minor axis, and 
position angle. The major and minor axes (FWHM) have 
been deconvolved using a 5.4" synthesized beam. Negative 
values indicate that the source size (before deconvolution) 
was smaller than the beam. These radio parameters arc 
described in further detail in WBHG. 

As demonstrated above, 95% of the counterparts within 
2" are actual FIRST identifications. There are, however, 
a number of reasons for studying the contents of the op- 
tical sky even farther from the FIRST positions: double 
and multiple-component sources are likely to have parent 
galaxies with substantial offsets from the radio positions, 
and the clustering of optical objects around FIRST sources 
appears to be a fruitful area for investigation (Helfand et 
al. 1998). In addition, catalog users may find it helpful to 
be able to calculate false rates for particular ranges of op- 
tical magnitude, radio flux density, or source morphology. 
Thus, we have constructed a second catalog (also avail- 
able on-line) which includes all optical objects appearing 
on either POSS-1 plate which lie within 20" of a FIRST 
position; 68% of the FIRST sources have at least one such 
counterpart. The format of this database is identical to 
that of the primary catalog; it contains nearly half a mil- 
lion entries representing 454,754 unique optical objects in 
the vicinity of 261,434 radio catalog entries. Also as the 
FIRST survey grows in size we shall periodically update 
these on-line optical catalogues. 

8. DISCUSSION 

The systematic optical identification of radio-selected 

samples with flux densities in the range 1-30 mJy began 
nearly twenty years ago using the first deep radio images 
obtained by Westerbork and the VLA. In Table 5, we sum- 
marize the published work on this problem. The area sur- 
veyed covers a little more than 50 deg^ and includes 1448 
radio sources between the limits of the all-sky, single dish 
surveys (~30 mJy) and the FIRST survey threshold (1.0 
mJy). A total of 455 sources have suggested identifica- 
tions, albeit to different levels of optical sensitivity, com- 
pleteness, and reliability. Only 270 counterpart candidates 
have been published for radio-selected objects at thresh- 
olds roughly equivalent to the POSS-I plates. The total 
number from this flux density range in the work reported 
here represents a sample of ~ 65,400 confidently identi- 
fied radio sources to this limit. Our project, then, repre- 
sents an increase of more than two orders of magnitude 
in the number of optically identified faint radio sources. 
What source classes are represented in this sample, and 
how do the results compare with our expectations from 
earlier work? 

The raw fraction of detections at the POSS-I limit, 
~ 18.3%, is similar to the mean value derived from previ- 
ous work (Table 5); the different optical bands used and 
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Fig. 18. — The fraction of FIRST sources with optical counterparts (matching to within 2") as a function of radio flux density. The shaded 
histogram includes objects classified as stars on at least one of the POSS-I plates; the remaining objects are classified as galaxies. 
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Fic;. 19. — The fraction of FIRST sources with stellar optical counterparts (matching to within 2") as a function of radio flux density. Given 
the tendency toward misclassification of faint galaxies as stellar, the true fraction of galaxy counterparts at faint flux levels is considerably 
higher than 60%. 
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Fig. 20. — Color-magnitude diagrams for the 7877 objects classified as stellar on both plates in the APM catalog, (a) Objects lying within 
1.5" of FIRST sources, (b) An equal number of stellar objects randomly selected from the same plates. The FIRST counterparts, fainter and 
bluer than the typical star, are mostly quasars. 
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Fig. 21. — Color-magnitude diagrams for the 16,047 objects classified as nonstellar on both plates in the APM catalog, (a) Objects lying 
within 4" of FIRST sources, (b) An equal number of nonstellax objects randomly selected from the same plates. The FIRST counterparts 
are, on average, brighter. 
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Fig. 22. — A contour and greyscale display of the fraction of objects classified as stellar on both plates which have radio counterparts. 
Darker regions have higher detected fractions. The fraction reaches nearly 10% for colors of O — E = 0. Contour levels are 0.0001, 0.0002, 
0.005, 0.001, 0.002, 0.005, 0.01, 0.02, and 0.05. 
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Fig. 23. — A contour and greyscale plot of the fraction of objects classified as nonstellar on both plates which have radio counterparts. 
More than 10% of all galaxies are detected at the brighter magnitudes. Contour levels are 0.002, 0.005, 0.01, 0.02, 0.05, 0.1, and 0.2. 
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the fuzzy lower bound to the POSS-I magnitudes easily 
explains the small discrepancy introduced when the cor- 
rections for multiple-component sources are applied (§5.4). 
The distribution amongst stellar and non-stellar counter- 
parts is also roughly similar; detailed comparisons arc not 
warranted given the inevitable difficulty of classifying ob- 
jects near the plate limits where most of the candidates 
lie. Further progress in quantifying the source populations 
represented, and in establishing their luminosity functions, 
requires CCD imaging and spectroscopy, work wc have be- 
gun to pursue for a variety of subsamples derived from this 
FIRST/ APM database. We discuss these briefly here. 

The comprehensive catalog of stellar radio observations 
by Wendkcr (1995) contains 18 stars in the FIRST suivey 
area detected at any radio frequency between 100 MHz 
and 30 GHz; of these, only a dozen have ever been seen 
at a 20 cm flux density in excess of 1 mJy. Our compre- 
hensive study of various stellar proper motion catalogs has 
produced an additional fourteen stellar radio source iden- 
tifications (Helfand et al. 1999). Proper motions remain 
the largest obstacle to our identification of stellar counter- 
parts fainter than Tycho catalog limit (~ 11*'* magnitude), 
and the advent of the GSC-II should significantly increase 
the number of stellar identifications to ~ 85. Nonethe- 
less, stars make a trivial contribution to the total radio 
source population above our 20 cm flux density threshold 
of 1 mJy. Other Galactic counterparts are similarly rare: 
we have detected only 3 (of the 8) known planetary neb- 
ulae and 7 (of the 18) known radio pulsars. In total, we 
expect that less than 0.03% of the FIRST catalog entries 
represent Galactic objects. 

There are, of course, a much larger number of opti- 
cal counterparts with stellar images on the POSS plates. 
These represent a mixture of quasars, BL Lac objects, 
AGN with sufficiently bright nuclei that they are clas- 
sified as stellar, H II galaxies with similarly bright nu- 
clei, and radio galaxies with sufficiently small diameters 
that they are unresolved. The FIRST Bright Quasar Sur- 
vey (White et al. 2000) has established the distribution 
amongst these categories for objects brighter than E=17.8: 
- 55% quasars, - 8% BL Lacs, ~ 9% AGN, - 4% radio 
galaxies, and ^ 17% H II galaxies, with the remaining ob- 
jects being a collection of stars and chance coincidences. 
We have begun a smaller survey to establish this distribu- 
tion at fainter magnitudes (Becker et al. 2001). 

The majority of the FIRST counterparts brighter than 
the POSS-I limit are galaxies. Two principal populations 
are represented: active galactic nuclei including both rela- 
tively nearby Seyferts and luminous elliptical hosts of more 
powerful radio galaxies up to z ~ 0.5, and a local popu- 
lation of star-forming galaxies (Condon 1992). Our work 
to date has consisted primarily of following up the coun- 
terparts to the several hundred bent, double- lobed radio 
galaxies which can act as tracers for galaxy clusters (Blan- 
ton et al. 2000, 2001). Studies of selected, magnitude- 
limited samples of the galaxy population are in progress. 

One important use of FIRST radio sources is as trac- 
ers of the general matter density of the Universe; their 
broad redshift distribution and high mean z samples a 
much larger volume than most surveys in other wave- 
length regimes. In applications such as determining the 
two-point (Cress et al. 1996), and higher order (Maglio- 



chctti et al. 1998) angular correlation functions, searching 
for a weak lensing signal on large angular scales (Refregier 
et al. 1999) and other such statistical studies, it may be 
useful to eliminate the low-redshift, star- forming galaxy 
population. The catalog presented here can be applied for 
such purposes. 

9. SUMMARY AND FUTURE WORK 

We have completed an optical identification program for 
FIRST radio sources by matching the FIRST catalog to 
a catalog constructed from the APM scans of the POSS-I 
survey. A detailed analysis of the properties of the APM 
catalog has been undertaken; we have recalibrated its pho- 
tometry using the APS, and have produced a factor of ~ 10 
improvement in astrometric accuracy by using the FIRST 
sources as a set of astrometric standards. We have dis- 
cussed in some detail the procedures for identifying radio 
source counterparts through positional coincidence, calcu- 
lating false match rates for various subsets of both cata- 
logs. For stellar counterparts of isolated radio sources, the 
matches are 90% reliable and 90% complete within an off- 
set of 1.2" ; over 98% of all matches within 1" are correct 
identifications. In total, nearly 73,000 radio source identi- 
fications arc recorded, an increase by a factor of ^ 200 over 
the total number of counterparts previously reported for 
radio-selected samples in the flux density range 1-30 mJy. 
We present color-magnitude diagrams for the radio coun- 
terparts and briefly discuss the various source populations 
represented. 

A catalog of all 454,754 optical objects within 20" of 
the 382,892 radio sources in the north Galactic cap re- 
gion FIRST has covered to date is available on the FIRST 
Webpage. As we complete radio observations and opti- 
cal matching for additional regions of the sky, this cata- 
log will be expanded accordingly. The results presented 
here suggest that future work to match FIRST with the 
GSC-II and the Sloan Digital Sky Survey (with which the 
FIRST sky coverage is matched) will yield enormous, valu- 
able samples of optically identified radio sources for use in 
a wide variety of applications from stellar radio astronomy 
to cosmology. 

The success of the FIRST survey is in large measure 
due to the generous support of a number of organizations. 
In particular, we acknowledge support from the NRAO, 
the NSF (grants AST-98-02791 and AST-98-02732), the 
Institute of Geophysics and Planetary Physics (operated 
under the auspices of the U. S. Department of Energy 
by Lawrence Livermore National Laboratory under con- 
tract No. W-7405-Eng-48), the Space Telescope Science 
Institute, NATO, the National Geographic Society (grant 
NGS No. 5393-094), Colmnbia University, and Sun Mi- 
crosystems. We thank Juan Cabanela for his assistance 
in matching with the APS catalog. RGM acknowledges 
support from the Royal Society. DJH is grateful for the 
support of the Raymond and Beverly Sackler Fund, and 
joins RHB and RLW in thanking the Institute of Astron- 
omy of the University of Cambridge for hospitality during 
much of this work. We thanks Mike Irwin for his explana- 
tory comments on a draft of this paper. This paper is 
Contribution Number 685 of the Columbia Astrophysics 
Laboratory. 



Optical IDs for 70,000 Radio Sources 



25 



REFERENCES 



Baadc, W., & Minkowski, R. 1954a, ApJ, 119, 206 
Baadc, W., & Minkowski, R. 1954b, ApJ, 119, 215 
Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559 
Becker, R. H., White, R. L.. Gregg, M., Laurent-Muehleisen, S., 
Brotherton, M., Impey, C, Chaffee, F., Richards, G. T., Heffand, 
D. J., Lacy, M., Courbin, F., & Proctor, D. D. 2001, ApJS, 135, 
227 

Benn, C. R., Grueff, G., Vigotti, M., & Wall, J. J. 1988, MNRAS, 
230, 1 

Blanton, E. L., Gregg, M. D., Helfand, D. J., Becker, R. H., & White, 

R. L. 2000, ApJ, 531, 118 
Blanton, E. L., Gregg, M. D., Helfand, D. J., Becker, R. H., & 

Leighly, K. M. 2001, AJ, 121, 2915 
Bolton, J. G., Stanley G. J., & Slee, O. B. 1949, Nature, 164, 101 
Bunclark, P. S., & Irwin, M. J. 1983, in Proc. Statistical Methods in 

Astronomy, ESA SP-201, edited by E. J. Rolfe, pl95 
Condon, J. J. 1992, ARAA, 30, 575 

Coleman, P. H., Condon, J. J., & Hazard, C. 1985, AJ, 90, 1437 
Condon, J. J., Condon, M. A., & Hazard, C. 1982, AJ. 87, 739 
Condon, J. J., & Mitchell, K. J. 1984, AJ, 89. 610 
Condon, J. J., Dickey J. M., & Salpeter, E. E. 1990, AJ 99, 1071 
Cress, C. M., Helfand, D. J., Becker, R. H., Gregg, M. D., & White, 

R. L. 1996, ApJ, 473, 7 
deRuiter, H. R., Zamorani, G., Parma, P., Hasinger, G., Hartner, G., 

Trumper, J., Burg, R., Giacconi, R., & Schmidt, M. 1997, A&A 

319, 7 

Evans, D. W. 1989, A&AS, 78, 249 

Fomalont, E. B., Windhorst, R. A., Kristian, J. A., & Kellerman, K. 

I. 1991, AJ, 102, 1258 
Georgakakis, A., Mobasher, B., Cram, L., Hopkins, A., Lidman, C. 

& Rowan-Robinson, M. 1999, MNRAS, 306, 708 
Gregg, M. D., Becker, R. H., White, R. L., Helfand, D. J., McMahon, 

R. G., & Hook, I. M. 1996, AJ, 112, 407 
Gunn, J. E. & Knapp, G. R. 1993, in Sky Surveys, ed. B. T. Soifer 

(San Francisco: ASP), 267 
Harrington, R. G. 1952, PASP, 64, 275 
Hayes, D. S., & Latham, D. W. 1975, ApJ, 197, 593 
Helfand, D. J., Yadigaroglu, L A., Berger, R., Postman, M., White, 

R. L., Lauer, T. R., Oegerle, W. R., & Becker, R. H. 1998, BAAS, 

30, 1407 

Helfand, D. J., Schnee, S., Becker, R. H., White, R. L., & McMahon, 

R. G. 1999, AJ, 117, 1568 
Helfand, D. J., Stone. R. P. S., Willman, B., White, R. L., Becker, 

R. H., Price, T., Gregg, M. D., & McMahon, R. G. 2001, AJ, 121, 

1872 

Hey, J. S. 1973, in "The Evolution of Radio Astronomy" (London: 

Elek Science), pl5 
Hjellming, R. M. 1998, in Galactic and Extragalactic Radio 

Astronomy, ed. G. L. Verschuur & K. L Kellerman (Berlin: 

Springer), p381 
Hog, E. et al. 1997, A&A 323, L57 

Hook, I. M, McMahon, R. G., Boyle, B. J. & Irwin, M. J. 1994, 

MNRAS, 268, 305. 
Hook, I. M., McMahon, R. G., Irwin, M. J., & Hazard, C. 1996, 

MNRAS, 282, 1274 
Hook, I. M., Becker, R. H., McMahon, R. G., & White, R. L. 1998, 

MNRAS, 297, 1115 
Hopkins, A. Mobasher, B., Cram, L., & Rowan-Robinson, M. 1998, 

MNRAS, 296, 839 
Humphreys, R. M., Landau, R. L., Ghigo, F. D., Zumach, W., 

LaBonte, A. E. 1991, AJ, 102, 395 
Irwin, M. J. 1985, MNRAS, 214, 575 
Irwin, M. 1994, lAU Commission 9 Newsletter, 5, 25. 



Johnston, K. et al. 1995, AJ, 110, 880 

Kron, R. G., Koo, D. C, & Windhorst, R. A. 1985, A&A, 146, 38 

Kibblewhite, E. J., Bridgeland, M. T., Bunclark, P. S. & Irwin, M. 
J. 1984, AstroiJoinicaJ Microdensitometry Conference NASA CP- 
2317, ed. Klinglcsmith, D. A. (Washington: NASA), p277 

Lu, N. Y., Hoffman, G. L., Salpeter, E. E., & Houck, J. R. 1996, 
ApJS, 103, 331 

Lund, J., & Dixon, R. 1973, PASP, 85, 230 

Maddox, S. J., Sutherland, W. J., Efstathiou, G. & Loveday, J., 

1990a, MNRAS, 243, 692 
Maddox, S. J., Efstathiou, G. & Sutherland, W. J. 1990b, MNRAS, 

246, 433. 

Magliocchetti, M., Maddox, S. J., Lahav, O., & Wall, J. V. 1998, 

MNRAS, 300, 257 
McMahon, R. G., & Irwin, M. J. 1992, in Digitised Optical Sky 

Surveys, eds. H. T. MacGillivray & E. B. Thomson (Dordrecht: 

Kluwer), p417 

Metcalfe, N., Fong, R., Shanks, T. 1995, MNRAS, 274, 769 
Minkowski, R. L., & Abell, G. O. 1963, "Stars and Stellar Systems 

Vol. HI, Basic Astronomical Data" (Chicago: University of 

Chicago Press), p481 
Mitchell, K. J., & Condon, J. J. 1985, AJ, 90, 1957 
Nilson, P. 1973, The Uppsala General Catalogue of Galaxies, 

(Uppsala Astronomical Observatory, Uppsala) 
Partridge, R. B., Hilldrup, K. C, & Ratner, M. I. 1986, ApJ, 308, 

46 

Patnaik, A. R., Browne, I. W. A., Wilkinson, P. N., Wrobel, J. M. 
1992, MNRAS, 254, 655 

Pennington, R. L., Humphreys, R. M., Odewahn, S. C, Zumach, W., 
& Thurmes, P. M. 1993, PASP, 105, 521 

Postman, M., Bucciajrelli, B. Sturch, C, Borgman, T. Casalegno, R., 
Doggett, J., & Costa, E. 1998a, in "New Horizons from Multi- 
Wavelength Sky Surveys (lAU Symposium 179), ed B. J. McLean, 
D. A. Golumbeck, J. J. E. Hayes, & H. E. Payne (Amsterdam: 
Kluwer), p379. 

Postman, M, Lauer, T. R., Szapudi, I., & Oegerle, W. 1998b, ApJ, 

506, 33 

Reid, I. N. et al. 1991, PASP, 103, 661 

Refregier, A., Brown, S. T., Kamionkowski, M. Helfand, D. J., 

Cress, C. M., Babul, A., Becker, R. H., & White, R. L. 1999, 

in "Cosmology with Wide Field Surveys", eds. S. Colombi & Y. 

Mellier (Paris: Editions Frontieres), p209 
Richards, E. A., Kellerman, K. I., Fomalont, E. B., Windhorst, R. 

A., & Partridge, R. B. 1998, AJ, 116, 1039 
Roscr, S., & Bastian, U. 1991, PPM Catalog (Heidelberg: Spektrum 

Akademisher Verlag) 
Ryle, M. Smith, F G., & Elsmore, B. 1950, MNRAS, 110, 508 
Snellen, I., McMahon, R. G., Hook, I. M., Browne, I. W. 2001, 

MNRAS, in press (preprint astro-ph/0110136) 
Urban, S. E., Corbin, T. E., Wycoff, G. L. 1998, AJ, 115, 2161 
Veron-Cetty, M.-P. & Veron, P. 1998, A Catalogue of Quasars 

and Active Nuclei (8th Edition), European Southern Observatory 

Scientific Report, 18, 1 
Wendker, H. J. 1995, A&AS, 109, 177 

Weistrop, D., Mall, J. V., Fomalont, E. B., & Kellerman, K. I. 1987, 
AJ, 93, 805 

White, R. L., Becker, R. H., Helfand, D. J., & Gregg, M. D. 1997, 

ApJ, 475, 479 
White, R. L. et al. 2000, ApJS, 126, 133 

Windhorst, R. A., Kron, R. G., & Koo, D. C. 1984, A&AS, 58, 39 
Windhorst, R. A., Miley G. K., Owen, F. N., Kron, R. G., & Koo, 
D. C. 1985, ApJ, 289, 494 



McMahon, White, Becker, & Helfand 



Table 1 

APM Background Source Density 



APM Source 


Popt 


2 


Properties 


(10-^ arcsec-2) 


(10-8 arcsec-^) 


All matches 


7.58 


7.10 


Stellar on both plates 


1.60 


0.06 


Stellar on at least 1 plate 


5.16 


3.85 


Non-stellar on both plates 


0.41 


0.13 


Non-stellar on at least 1 plate 


2.98 


3.57 


Detected on both plates 


2.56 


0.51 


Detected on only 1 plate 


5.02 


6.59 


Detected only on red plate 


2.60 


2.76 


Detected only on blue plate 


2.43 


3.83 



Table 2 

Confidence Radii and False Rates for Various Source Types 



APM Sources 


FIRST Sources 


90% Radius'' 


Number of 


Number of Chance 


Percentage of 






(arcsec) 


True Matches^ 


Coincidences'^ 


False Matches 


All 


All 


4.2 


65400 


13000 


16.6% 


Non-stellar** 


All 


4.0 


22600 


640 


2.7% 


Stellar** 


All 


1.7 


7600 


520 


6.3% 


All 


Unresolved 


2.0 


17800 


850 


4.6% 


Non-stellar <* 


Unresolved 


2.2 


6200 


60 


0.9% 


Stellar** 


Unresolved 


1.1 


2300 


60 


2.4% 



'^APM-FIRST separation that encloses 90% of the true associations. 

''Estimated number of true matches within the 90% radius. 

*^Estimated number of chance background matches within the 90% radius. 

•^APM classification on both E and O plates. 



Table 3 

Confidence and Reliability for Single Band Detections 





Blue-only objects 


Red-only objects 


Radius 


Completeness 


Reliability 


Completeness 


Reliability 


(arcsec) 










1.0 


0.56 


0.93 


0.46 


0.98 


2.0 


0.78 


0.82 


0.70 


0.94 


3.0 


0.88 


0.71 


0.78 


0.89 


4.0 


0.94 


0.60 


0.84 


0.83 


5.0 


0.97 


0.50 


0.87 


0.77 
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Table 4 

Sami'LK Pace fkom FIRST-APM Soruci; C'atalcx! 



RA 




Dec 


ARA 


ADec 


Sep 


Eel 


Epsf 


E 


Ocl 


Opsf 


o 


O- 


E 


Plate 


(14) 


S, 


Maj 


Mill 


PA 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(15) 


(16) 


(17) 


(18) 


08 00 00.099 


+32 


06 09.80 


-1.86 


-0.23 


1.87 


1 


2.49 


20.72 


-1 


-0.31 


21.88 


1 


.16 


1344 


1.10 


1.09 


3 


.0 


-2.7 


173.9 


08 00 00.234 


+43 


45 30.57 


-0.56 


-0.21 


0.59 


-1 


-0.07 


20.22 





0.00 


-22.17 


-9, 


.99 


1329 


8.52 


9.27 


2 


.1 


1.0 


105.7 


08 00 00.243 


+22 


26 24.80 


0.56 


0.38 


0.67 


1 


5.09 


19.77 





0.00 


-21.85 


-9 


.99 


226 


1.21 


1.05 


1, 


.5 


-3.0 


179.1 


08 00 00.397 


+44 


39 14.16 


2.00 


0.17 


2.00 





0.00 


-20.76 


-1 


0.24 


21.31 


-9 


.99 


1317 


19.48 


22.0 


2, 


.3 


1.5 


80.0 


08 00 00.752 


+32 


46 50.87 


0.14 


-0.77 


0.78 


1 


4.60 


19.54 


2 


3.11 


22.46 


2 


.91 


989 


4.20 


3.90 


-1 


.1 


-1.7 


67.1 


08 00 01.358 


+47 


01 38.37 


-0.16 


-0.16 


0.23 





0.00 


-20.76 


-1 


0.97 


22.23 


-9 


.99 


1317 


2.68 


3.12 


2 


.2 


2.2 


90.3 


08 00 02.263 


+51 


36 09.10 


-0.31 


0.30 


0.43 


1 


13.85 


12.08 


2 


47.05 


15.00 


2, 


.92 


985 


1.49 


1.75 


4, 


.2 


-2.0 


36.6 


08 00 02.302 


+31 


17 20.98 


2.27 


1.38 


2.65 


-1 


-0.11 


20.48 


1 


2.85 


22.21 


1 


.73 


1344 


2.16 


2.73 


3 


.3 


2.1 


169.9 


08 00 02.515 


+35 


49 27.12 


-2.02 


4.07 


4.54 


-1 


-1.52 


18.88 


-1 


-0.96 


21.40 


2, 


.52 


989 


1.09 


0.87 


2, 


.2 


-3.6 


132.4 


08 00 02.702 


+22 


48 54.11 


-0.79 


-0.66 


1.03 


2 


24.52 


17.29 


1 


3.15 


21.62 


4 


.33 


226 


1.87 


2.91 


5, 


,3 


2.6 


94.6 


08 00 02.932 


+27 


17 27.51 


-0.64 


0.46 


0.79 


-1 


0.80 


20.48 


-1 


0.40 


20.77 





.29 


1344 


5.37 


5.56 


1 


.3 


0.6 


0.4 


08 00 03.250 


+36 


29 42.11 


0.46 


0.14 


0.48 


1 


3.68 


19.67 





0.00 


-22.18 


-9, 


.99 


989 


46.71 


58.72 


3, 


.7 


1.5 


21.6 


08 00 03.460 


+50 


30 05.76 


0.39 


-0.16 


0.42 


-1 


2.15 


20.23 


1 


2.61 


22.41 


2 


.17 


985 


1.07 


0.98 





.5 


-2.2 


0.8 


08 00 04.050 


+23 


26 16.30 


-2.87 


1.68 


3.32 


2 


19.81 


9.92 


2 


21.84 


11.49 


1, 


.56 


226 


6.00 


6.23 


1 


.1 


1.0 


180.0 


08 00 04.378 


+46 


09 59.44 


0.83 


0.08 


0.84 


1 


7.16 


18.45 


-1 


1.89 


22.31 


3 


.86 


1317 


2.01 


2.76 


4, 


.2 


2.3 


134.0 


08 00 04.981 


+28 


06 08.29 


-2.43 


0.54 


2.49 


2 


67.99 


14.50 


1 


23.40 


18.21 


3 


.71 


1344 


1.02 


1.55 


4, 


.1 


3.8 


90.5 


08 00 05.237 


+50 


10 23.62 


-0.44 


0.16 


0.47 


-1 


1.15 


18.93 





0.00 


-22.28 


-9, 


.99 


1317 


3.09 


3.88 


5, 


.3 


-2.4 


60.1 


08 00 05.463 


+32 


11 23.10 


0.43 


0.77 


0.88 


1 


7.23 


20.30 





0.00 


-21.83 


-9 


.99 


1344 


2.92 


3.14 


2, 


.5 


-1.2 


23.0 


08 00 06.509 


+43 


09 53.78 


-1.72 


2.62 


3.14 


2 


8.95 


20.09 





0.00 


-22.17 


-9 


.99 


1329 


2.69 


2.07 


-2, 


.3 


-2.8 


138.3 


08 00 06.555 


+32 


01 47.98 


0.85 


0.40 


0.94 


2 


10.27 


20.44 





0.00 


-21.83 


-9 


.99 


1344 


1.12 


0.72 


-2 


.0 


-3.9 


144.3 


08 00 07.243 


+25 


22 19.95 


0.36 


-0.19 


0.41 


-1 


0.29 


19.99 


-1 


0.09 


20.40 





.40 


226 


20.86 


38.93 


7, 


.7 


2.0 


127.2 


08 00 09.613 


+38 


06 50.00 


0.52 


0.72 


0.89 


1 


2.07 


21.03 





0.00 


-22.18 


-9 


.99 


989 


1.05 


1.39 


5 


.5 


-2.0 


141.9 


08 00 09.892 


+44 


36 45.86 


-1.85 


2.50 


3.11 


-1 


0.11 


20.60 





0.00 


-22.28 


-9 


.99 


1317 


1.21 


1.36 


6, 


.5 


-3.7 


24.4 


08 00 10.262 


+41 


20 36.84 


-1.42 


-1.72 


2.23 





0.00 


-20.84 


1 


8.15 


22.12 


-9 


.99 


1329 


1.36 


0.80 


-3 


.3 


-3.6 


91.9 


08 00 10.571 


+34 


02 20.90 


-0.19 


-0.93 


0.95 


-1 


1.24 


19.27 


2 


6.19 


22.30 


3 


.03 


989 


1.60 


3.45 


8, 


.9 


2.7 


180.0 


08 00 10.974 


+24 


09 03.42 


-3.57 


-2.23 


4.21 


1 


7.99 


18.30 


1 


7.09 


20.28 


1 


.98 


226 


2.31 


2.34 


2, 


.5 


-2.1 


6.0 


08 00 11.108 


+33 


14 27.82 


0.03 


-0.71 


0.71 


1 


2.39 


21.05 





0.00 


-22.18 


-9 


.99 


989 


11.21 


11.67 


1 


.4 


0.7 


4.3 


08 00 11.197 


+30 


31 48.00 


-0.28 


-1.17 


1.21 


1 


3.24 


20.41 





0.00 


-21.83 


-9 


.99 


1344 


26.07 


25.93 





.4 


-0.7 


110.6 


08 00 12.213 


+22 


49 27.46 


-0.02 


1.88 


1.88 


-1 


0.29 


19.97 





0.00 


-21.85 


-9 


.99 


226 


1.37 


2.14 


4, 


.3 


3.9 


87.2 


08 00 12.515 


+26 


22 16.93 


-2.49 


-0.59 


2.56 





0.00 


-20.26 


1 


2.03 


20.53 


-9 


.99 


226 


26.25 


27.72 


1 


.3 


1.3 


180.0 


08 00 12.603 


+45 


39 23.51 


3.34 


-3.24 


4.66 


2 


7.07 


20.19 





0.00 


-22.28 


-9 


.99 


1317 


1.97 


3.15 


4, 


.3 


4.1 


0.2 
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Table 5 

Previous Optical Identifications of Faint Radio Sources 



Position 


Area 


5(1.4 GHz) 


Sources 


Suggested Identificatio] 


as 


Reference'' 


(19.")0) 


((l0g2) 


(u).Ty) 


(1 -M) m.T\-) 


POSS-1 


riijiit w 


Total 




1204+1130 


9.0 




102 


21 


- 


21 


1 


130D+2900 
















1 71 s-i-i^nnn 

1 ( IOt^OUUU 


o.o 


u.o 




45 


100 




9 


nm A 1 1 ^Qfi 
















Uo4U+444U 
















0852+1716 


0.44^^ 


0.1 


37 


11 





11 


3 


1300+3034 


0.44 


0.1 


40 


7 


5 


18 


4 




b 


0.6 


66 


11 


14 


25 


5 


0854+1730 


9.0 


6.0 


130 


31 





31 


6 


1204+1130 
















0310+8008 


0.13 


1.0= 


14 


3 





3 


7 


1635+6620 


0.13 














0015+1533 


0.01 


0.3= 


13 


1 


3 


4 


8 


NGP 


12.57 




168 


31 


35 


66 


9 


0927+4710 


12.0 




329 


53 


4 


57 


10, 11 


1416+5242 


0.5 


0.03= 


6 


1 


1 


2 


12 


HDF 


0.02 


0.02^1 


5 





4 


4 


13 


Lockman Hole 


0.35 


0.12 


28 


5 


1 


6 


14 


Phoenix Deep 


3.0 


0.10 


270 


50 


12 


62 


15, 16 


Totals: 


53.1 




1448 


270 = 18.6% 




455 




This work: 


4150 


1.0 


357.985 


65.400'= 18.3% 




65,400 





''Partial overlap with the sample of reference 6 
''Partial overlap with sample of reference 2 
=6 em sm'vey; fluxes translated to 20 em using a < —0.5 
■^^'3.6 cm survey; fluxes translated to 20 cm using a < —0.5 

"References: 1. Condon, Condon, & Hazard 1982; 2. Windhorst, Kron, & Koo 1984; Kron, Koo, and 
Windhorst 1985; 3. Condon & Mitchell 1984; 4. Mitchell & Condon 1985; 5. Windhorst et al. 1985; 6. 
Coleman, Condon, & Hazard 1985; 7. Partridge, Hilldrup & Ratner 1986; 8. Weistrop et al. 1987; 9. Benn et 
al. 1988; 10. Condon, Dickey, & Salpeter 1990; 11. Lu et al. 1996; 12. Fomalont et al. 1991; 13. Richards et 
al. 1998; 14. deRuiter et al. 1997; 15. Hopkins et al. 1998; 16. Georgakakis et al. 1998. 

^65,400 is the number of 90%-confident identifications given in Table 2; it does not include lower confidence 
associations or objects identified with double radio sources, which boost the number of suggested identifications 
to ~ 73, 000. 



